Many-body aspects of coherent atom-molecule oscillations 
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We study the many-body effects on coherent atom-molecule oscillations by means of an effective 
quantum field theory that describes Feshbach-resonant interactions in Bose gases in terms of an 
atom-molecule hamiltonian. We determine numerically the many-body corrections to the oscillation 
frequency for various densities of the atomic condensate. We also derive an analytic expression that 
approximately describes both the density and magnetic-field dependence of this frequency near the 
resonance. We find excellent agreement with experiment. 
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Introduction — One of the most remarkable applica- 
tions of Feshbach resonances in doubly spin-polarized al- 
kali gases P, 0) 01 is the observation of coherent atom- 
molecule oscillations In this last experiment Donley 
et al. used the Feshbach resonance at Bq ~ 154.9 Gauss 
in the |/ = 2;m/ = —2) state of ^^Rb to perform a 
Ramsey- type experiment, consisting of two short pulses 
in the magnetic field towards resonance separated by a 
longer evolution time. As a function of this evolution 
time an oscillation in the number of condensate atoms 
was observed. Over the investigated range of magnetic 
field during the evolution time, the frequency of this oscil- 
lation agreed exactly with the molecular binding energy 
found from a two-body coupled-channels calculation 
indicating coherence between atoms and molecules. 

Very recently, Claussen et al. have performed a sim- 
ilar series of measurements over a larger range of mag- 
netic fields j^. It was found that close to resonance the 
frequency of the oscillation deviates from the two-body 
molecular binding energy, which indicates that many- 
body effects play an important role in this regime. It 
is the main purpose of this Letter to present the theory 
that explains this deviation quantitatively. 

The first mean-field theory for Feshbach-resonant in- 
teractions in Bose-Einstein condensed gases is due to 
Drummond et al. and Timmermans et al. jsj and 
introduces the physical picture of an interacting atomic 
condensate coupled to a noninteracting molecular con- 
densate. Although this theory contains the correct reso- 
nant scattering amplitude for the atoms, it does not con- 
tain the correct molecular binding energy. By studying 
the fluctuations around this mean-field theory within the 
Hartree-Fock-Bogoliubov approximation, it is possible to 
also incorporate the correct binding energy |^, llfll |. 
This comes about because the latter approach explicitly 
contains a so-called anomalous density, or pairing field, 
for the atoms which, after elimination, leads to a shift 
in the coupling constants and the molecular binding en- 
ergy. Unfortunately, however, the elimination does not 
lead to a proper renormalization of all the coupling con- 
stants. In particular the interaction between condensate 



atoms and non-condensate atoms, which is ultimately re- 
sponsible for the many-body corrections to the molecu- 
lar binding energy that are of interest to us here, is not 
correctly described within the Hartree-Fock-Bogoliubov 
approximation. 

Relying on the anomalous density for the description of 
the molecular properties also makes the theory inapplica- 
ble above the critical temperature for Bose-Einstein con- 
densation, since the anomalous density is proportional 
to the atomic condensate density. To overcome all these 
problems, it is convenient to formulate an effective quan- 
tum field theory that incorporates the exact two-body 
physics not at the mean-field level but at the quantum 
level and is hence applicable both above and below the 
critical temperature. We have recently derived such an 
effective quantum field theory by starting from the mi- 
croscopic hamiltonian for an atomic gas with a Feshbach 
resonance and explicitly summing all the ladder diagrams 
[llj |. Here, we apply the mean-field theory for the Bose- 
Einstein condensed phase of the gas that results from this 
effective quantum field theory, to the study of the recent 
experiments by Claussen et al. |^. 

Atom-molecule coherence — The mean-field equations 
consist of coupled equations for the macroscopic wave 
function ?/;a(x, i), describing the atomic Bose-Einstein 
condensate, and the macroscopic wave function '0m(x, t) 
that describes the condensate of bare molecules, i.e., with 
the appropriate bound-state wave function in the closed 
channel potential of the Feshbach problem. The mean- 
field equations for the coupled atom-molecule system are 
ultimately given by 
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where T^^ = 47rabgS.^/m is the off-resonant two-body 
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T(ransition) matrix with abg the off- resonant background 
scattering length and m the mass of one atom. The atom- 
molecule coupling g is found from experiment by adiabat- 
ically eliminating the molecular wave function and using 
the fact that the resulting magnetic-field dependent scat- 
tering length of the atoms must be equal to 



Close to the resonance it thus follows that 
em(S) — —h'^/[ma{B)'^], which is the correct molecular 
binding energy in vacuum Note that in the absence 
of the molecular self energy the binding energy would 
be equal to the detuning, which is incorrect close to 
resonance. Moreover, the residue of the pole is given by 



a{B) = flbg 1 
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with AB and Bq the experimental width and position 
of the resonance, respectively. This procedure results in 
g = hy/2TTahgABA^/m, where we have made use of the 
fact that the detuning of the bare molecular state is given 
by S{B) = A^(B - Bq) with A^ ~ -2.2^b for ^^Rb 5] 
and fiB the Bohr magneton. 

The mean-field equation for the molecular condensate 
contains a fractional derivative, corresponding to the re- 
tarded self energy of the molecules. In momentum and 
frequency space, this self energy reads hY^[^^ (k, lu) — 

-{g^m^/^/2Trh^)i^hLU - h^k^/4m - 2nSHF U\. The 
square-root behavior is a result of the Wigner threshold 
law for the decay of a molecule with total energy huj and 
center-of-mass momentum Kk into the two-atom contin- 
uum. Due to the mean-field interaction of the noncon- 
densed atoms with the condensate, the decaying molecule 
has to overcome a mean- field barrier given by 2hY]^^ . 
Here, denotes the Hartree-Fock self energy for the 

noncondcnscd atoms. Neglecting the momentum depen- 
dence of this self energy, it is given by 
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where ria = iV'aP is the density of the atomic conden- 
sate, fi its chemical potential, and we have used that a 
collision between a condensed atom and a noncondensed 
atom has a mean-field shift oi 11 + Far from res- 

onance the energy-dependence of the interactions can be 
safely ignored and the Hartree-Fock self energy becomes 
equal to 8'Ka{B)h^ns./m, as expected. 

The molecular binding energy is given by the pole of 
the molecular propagator at zero momentum, which from 
Eq. is seen to be given by 

GW(0,c.) 



hLj + iO-5{B) + 



(4) 

In the limit of vanishing condensate density Ug, and for 
negative detuning, it has a pole at 
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and always smaller than one. Physically, the latter can 
be understood from the fact that the dressed molecular 
bound state near the Feshbach resonance is given by 



I Xm; dressed) 
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-k;open) , (7) 



where the coefficient C(k) denotes the amplitude of the 
dressed molecular state to be in the open channel of 
the Feshbach problem and the two atoms having mo- 
menta k and — k, respectively. They are normalized as 
/ dk|C(k)|V(27r)3 = 1 - Z{B). The dressed molecular 
state therefore only contains with an amplitude \J Z{B) 
the bare molecular state |xm;bare). Close to resonance 
we have that Z{B) <C 1, whereas it approaches one far off 
resonance. With respect to this remark it is important to 
note that the result of the Hartree-Fock-Bogoliubov the- 
ory for the density of the molecular condensate should be 
multipHed by a factor 1/Z{B) :s> 1 to obtain the density 
of real dressed molecules, since in this theory always the 
density of bare molecules is calculated 0, ^ . 

Collective modes — To study the many-body effects on 
the frequency of the coherent atom-molecule oscillations 
it is important to realize that these oscillations are in 
fact a collective mode where the atomic condensate den- 
sity oscillates out-of-phase with the molecular conden- 
sate density. It is thus worthwhile to study the collective 
modes of the mean-field equations in Eq. (Q) and look for 
solutions of the form 



(8) 



After substitution into the mean-field equations in 
Eq. Q the eigenmodes are found by diagonalizing the 
resulting 4x4 matrix. This yields a dispersion rela- 
tion with two branches, one corresponding to the gap- 
less Bogoliubov modes and one that corresponds to the 
atom-molecule oscillations. The zero-momentum part of 
the latter corresponds to the experimentally observed fre- 
quency of the coherent atom-molecule oscillations. Note 
that for this calculation the evaluation of the molecular 
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FIG. 1; Dispersion relation for the collective modes of the 
atom-molecule system for an atomic condensate density of 
na = 2 X 10^^ cm"^, at a magnetic field of B = 156 G. The 
momentum is indicated in units of 1/^, where ^ — l/y/T6nan2 
is the coherence length. The upper solid line shows the gap- 
less dispersion relation for phonon-like excitations. The lower 
solid line indicates the dispersion for atom-molecule oscilla- 
tions. The dashed line shows the effective Bogoliubov disper- 



self energy can be performed exactly because every term 
in the right-hand side of Eq. is an eigenfunction of 
the operator under the square root. 

In Fig. n we show the dispersion relations found by 
means of the above procedure, for a fixed atomic conden- 
sate density of = 2 x 10^^ cm~'^ at a fixed magnetic 
field of _B = 156 G. Physically, the upper branch corre- 
sponds to the phonon-likc excitations in the atomic con- 
densate and the lower branch to coherent atom-molecule 
oscillations. The dashed line denotes the Bogoliubov dis- 
persion for the scattering length a{B). At low momenta 
the phonon branch corresponds with the Bogoliubov dis- 
persion, as expected. At higher momenta the dispersion 
starts to deviate from the Bogoliubov result, due to the 
energy and momentum dependence of the resonant inter- 
actions that reduce the scattering amplitude. The disper- 
sion of the lower branch obeys fiuj\^ ~ —Tt-loj + h'^'k'^ /Am, 
where hujj is the Josephson frequency that is observed 
by Claussen et al. in their Ramsey experiment. This 
dispersion is negative due to the fact that we are dealing 
with a metastable situation. For negative detuning the 
true ground state contains almost all atoms in the form 
of molecules. 

In Fig. 121 the results are shown for the frequency of 
the atom-molecule oscillations as a function of the mag- 
netic field, for three different densities. Clearly, for in- 
creasing density the frequency starts to deviate from the 
two-body result. In addition, we show in Fig.|21the fre- 
quency of the atom-molecule oscillations relative to the 
two-body binding energy, as a function of the atomic con- 
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FIG. 2: Frequency of coherent atom-molecule oscillations. 
The solid lines show the frequency as a function of the mag- 
netic field for three different densities of the atomic conden- 
sate. The dashed line shows the two-body result for the molec- 
ular binding energy. The experimental points are taken from 
Ref. i^. 



densate density. The calculation is performed in this case 
for several values of the magnetic field. Both the mag- 
netic field and the atomic density dependence of the fre- 
quency can be understood as follows. We first observe 
that for the range of magnetic fields and atomic densi- 
ties that are explored experimentally, the difference be- 
tween the energy of the dressed molecular state and the 
threshold of the two-particle continuum is almost inde- 
pendent of density and equal to eni(S). It is interesting 
to note that only this difference is independent of density, 
whereas both quantities individu ally show a substantial 
mean- field shift of about 2hJ:^'^ m|. As a result of the 
cancellation, however, also the wave-function renormal- 
ization factor Z is almost density independent and equal 
to Z{B). Expressing our mean-field equations in Eq. 
in terms of the condensate wave function for the dressed 
molecules by replacing Vm(x,i) by ^/Z{B)ipin{ji,t), we 
see that the coupling between the atomic condensate and 
the dressed molecular condensate is reduced by a factor 
of Z{B) <C 1 0. We thus expect the frequency to 
obey 



(9) 



This analytic expression indeed turns out to give a first 
approximation to the deviation of the frequency from the 
two-body result. 

In order to confront our results with the experimen- 
tal data we have to realize that the experiments are per- 
formed in a magnetic trap. Taking only the ground-states 
into account for both the atomic and the molecular con- 
densates, this implies effectively that the atom-molecule 
coupling g is reduced by an overlap integral. Hence we 
define the effective homogeneous condensate density by 
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FIG. 3: Frequency of coherent atom-molecule oscillations as 
a function of the atomic condensate density, relative to the 
two-body binding energy. The solid lines show the result for 
different magnetic fields. 



means of ria = A^a [/ '^x0^(x)0m(x)] ^, where TVa denotes 
the number of condensed atoms and (/'a(x) and (j)m{^) 
denote the atomic and molecular ground state wave func- 
tion, respectively. For the experiments of Claussen et al. 
this results in an effective density of ria — 2 x 10^^ cm~^ 
1^. Fig. 121 clearly shows an excellent agreement with the 
experimentally observed frequency for this density. 

It is important to note that there are two hidden as- 
sumptions in the above comparison. First, we have used 
that the dressed molecules are trapped in the same exter- 
nal potential as the atoms. This is not obvious because 
the bare molecular state involved in the Feshbach res- 
onance is high-field seeking and therefore not trapped. 
However, Eq. © shows that near resonance almost all 
the amplitude of the dressed molecule is in the low-field 
seeking open channel and its magnetic moment is there- 
fore almost equal to twice the atomic magnetic moment. 
Second, we have determined the frequency of the coher- 
ent atom-molecule oscillations in equilibrium. In con- 
trast, the observed oscillations in the number of conden- 
sate atoms is clearly a nonequilibrium phenomenon. This 
is, however, expected not to play an important role be- 
cause the Ramsey-pulse sequence is performed on such a 
fast time scale that the response of the condensate wave 
function can be neglected. 

Conclusions — With the linear-response calculation 
presented in this Letter we have obtained excellent agree- 
ment with the experimental results on the frequency of 
coherent atom-molecule oscillations. The next step is a 
more detailed understanding of other quantities that are 
of interest for the two-pulse experiments . This in- 
cludes a quantitative study of the number of condensed 
atoms as a function of time that goes beyond the linear 
approximation discussed here. Experimentally an over- 



all decay of the number of condensed atoms is observed 
and also that the oscillations have a finite damping rate. 
Both effects increase as one approaches the resonance. 
We expect that an important contribution to these ef- 
fects is due to the so-called rogue-dissociation process , 
which has its physical origin in the decay of a molecule 
into two noncondensed atoms that in the center-of-mass 
system have opposite momenta k and — k, respectively. 
In equilibrium this process is always forbidden due to 
energy conservation, since the molecular state lies be- 
low the two-atom continuum threshold. However, in the 
nonequilibrium setting of the experiments of interest, it 
may occur due to the fact that the detuning is strongly 
time-dependent. In the mean- field equations in Eq. it 
is the fractional derivative term corresponding to the self 
energy of the molecules that automatically incorporates 
the rogue-dissociation effect. A full numerical simulation 
of these equation is challenging due to the fact that we 
are essentially dealing with a term that is nonlocal in 
time. Nevertheless, work in this direction is in progress 
and will be reported in a future publication. 
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The full expression for the self energy of the molecules 
is given by REj^^(k,tj) = -(g^m^^'^ /2nn^)i^{l + 
iahgy^/ny^, with z ^huj- 7i^k^/4m - 2?iE™. Close 
to resonance only the square-root in the numerator is im- 
portant. We will use this approximation throughout the 
paper. 



